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; ABSTRACT 
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^ We present a study of correlations between X-ray spectral and timing prop- 

' erties observed from a number of Galactic Black Hole (BH) binaries during hard- 

K ■ soft state spectral evolution. We analyze 17 transition episodes from 8 BH sources 

■ observed with Rossi X-ray Timing Explorer (RXTE). Our scaling technique for 
)> I' 

I ' BH mass determination uses a correlation between spectral index and quasi- 

^ ■ periodic oscillation (QPO) frequency. In addition, we use a correlation between 

^ ■ index and the normalization of the disk "seed" component to cross-check the BH 

mass determination and estimate the distance to the source. While the index- 

T— I ■ QPO correlations for two given sources contain information on the ratio of the BH 

^ ■ masses in those sources, the index-normalization correlations depend on the ratio 

•/^ ■ of the BH masses and the distance square ratio. In fact, the index-normalization 
00 ■ 

^ ■ correlation also discloses the index-mass accretion rate saturation effect given 

^ ■ that the normalization of disk "seed" photon supply is proportional to the disk 

^ . mass accretion rate. We present arguments that this observationally established 

O ■ index saturation effect is a signature of the bulk motion (converging) flow onto 

. black hole which was early predicted by the dynamical Comptonization theory. 

^ . We use GRO J1655-40 as a primary reference source for which the BH mass, 

. distance and inclination angle are evaluated by dynamical measurements with 
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unprecedented precision among other Galactic BH sources. We apply our scal- 
ing technique to determine BH masses and distances for Cygnus X-1, GX 339-4, 
4U 1543-47, XTE J1550-564, XTE J1650-500, H 1743-322 and XTE J1859-226. 
Good agreement of our results for sources with known values of BH masses and 
distance provides an independent verification for our scaling technique. 

Subject headings: accretion, accretion disks — black hole physics — stars:individual 
(Cygnus X-1), individual (GRO J1655-40), individual (XTE J1550-564), individ- 
ual (H 1743-322), individual (GX 339-4), individual (XTE 1650-500), individual 
(XTE 1859+226), individual (4U 1543-47) :radiation mechanisms: nonthermal — 
physical data and processes 



1. Introduction 



Determination of masses of Galactic black holes (BHs) is one of the most important 
tasks in modern astronomy. In general, the knowledge of the BH mass distribution of the 
Galaxy can provide important clues on stellar evolution. It also can con strain the maximum 



mass of a neutron star and the minimum mass of a BH. For example, iRhoades fc Ruffini 



(I1974J ) predicted that the maximum mass of neutron star, formed as a result of a supernova 
explosion, is probably about 3.2 solar masses. Study of the characteristic changes in spectral 
and variability properties of X-ray binaries proved to be a valuable source of information on 
the physics governing the accretion processes and on the fundamental parameters of BHs. 

BH obs ervational appearance is conventionally descr i bed in terms of BH state classi- 



fication (see Remillard fc McChntock 2006; Belloni 2005; Klein- Wolt &: van der Klis 2008 



for different flavors of BH states definitions). We adopt the following general BH state clas- 
sification for five major BH states: gmescent, low-hard (LHS), intermediate (IS), high-soft 
(HSS) and very high states (VHS). When a BH transient goes into outburst it leaves a qui- 
escent state and enters LHS, a low luminosity state with the energy spectrum dominated 
by a Comptonization component combined (convolved) with a weak thermal component. 
The photon spectrum in LHS is presumably a result of Comptonization (upscattering) of 
soft photons, origin ated in a relatively weak accr etion disk, off electrons of the hot ambi- 
ent plasma [see e.g. ISunyaev fc Titarchuk I (jl980l ))]. Variability in LHS is high (fractional 
root-mean-square variability is up to 40%) and presented by a fiat-top broken power law 
(white-red noise) shape, accompanied by quasi-periodic oscillations (QPOs) in the range 
of 0.01-30 Hz, observed as narrow peaks in the power density spectrum (PDS). The HSS 
photon spectrum is characterized by a prominent thermal component which is probably a 
signature of a strong emission coming from a geometrically thin accretion disk. A weak 



- 3 - 



power-law component is also present at the level of not more than 20% of the total source 
flux. In the HSS the fiat-top variability ceases, QPOs disappear and PDS acquires a pure 
power-law shape. The total variability in HSS is usually about 5% fractional rms. The IS is 
a transitional stage between LHS and HSS. 

Timing and spectral properties of an accreting BH are tightly correlated (see a com- 
prehensive set of references in Table [1]) . Correlations between the spectral hardness (photon 
index) and the characteristic frequencies of the quasi-periodic variability observed in the 
l ightcurves of BH sources h as been proposed f or use as a tool to determin e a B H mass 



(ITitarchuk. fc Fiorito 1 12004| . hereafter TF04). IShaposhnikov fc Titarchuk I (120071 ). here- 
after ST07, employed this idea to measure the mass of a prominent BH source Cygnus 
X-1 and obtained a BH mass of Mcygxi = (8.7 ± O.8)M0, which has a better preci- 
sion than conventional dynamical methods. Scaling argum ents were also used toesti- 



Fiorito fc Titarchuk 



mate masses in a number of ult r a-lurn i nous X-ray sources 

Dewangan. Titarchuk fc Griffiths I (120061 ): IStrohmayer et al. I (j2007l ). In this Paper we fur- 
ther explore this technique by analyzing more representative data sets from Galactic BH 
sources. We concentrate our efforts on the study of correlations between spectral index 
and characteristic sub-second variability frequencies and the total accretion disk luminos- 
ity. These two correlation patterns carry the most direct information on the BH mass 
and source distance. We also investigate the possibility that the shape of the correlation 
pattern can provide a direct signature of the bulk moti on (converging) flow onto a black 
hole, or, in other words, the signature of a black hole [see ITitarchuk fc Zannias I (119981 ) and 
Laurent fc Titarchuk I (119991 ) for more details on this subject]. We present observational and 
theoretical arguments that the index-mass accretion rate saturation effect is a signature of 
a converging flow, when the mass accretion rate exceeds the Eddigton limit, which can only 
exist in BH sources. 

In the presented study we enhance the scaling method by considering joint scaling of 
spectral index both in QPO frequency and disk component normalization (disk flux) do- 
mains. In the LHS and IS only a small part of the disk emission component is seen directly. 
The energy spectrum is dominated by a Comptonized component seen as a power law. To 
calculate the total normalization of the "seed" disk black body component we model the 
spectrum with a Generic Comptonization model which consistently convolves a disk black 
body with the Green's (response) function of the Compton Corona to produce the Comp- 
tonized component (see details in §2). Calculated in this way the spectrum normalization 
is directly related to the BH mass, emission efficiency, source distance and geometry. This 
allows us to apply the total disk normalization as an additional scaling variable and to 
impose an additional constraint on the BH mass and distance ratios. Thus, we use the 
index-QPO and index-disk normalization correlations to constrain the BH masses and dis- 
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tances for Galactic BH sources H 1743-322, GX 339-4, XTE J1550-564, XTE J1650-500, 4U 
1543-47, XTE J1859+226 and Cygnus X-1 based on their expected dependences on the mass 
of the central BH, the source distance and geometry. As our reference values for the scaling 
measurements, we use the previously measured BH mass, system distance and inclination 
for Galactic micro-quasar GRO J1655-40. 

The description of RXTE data for each source is given in ^ We provide the details of 
state transitions analysis in §3] and we present the results of the scaling method applied to 
the data in §U We present a physical picture of the accretion process which is consistent 
with the observed phenomenology of state transitions in ^ Specifically, it provides a theo- 
retical basis and an explanation for the observed scaling index-QPO and index-normalization 
pattern s, in Appenix Rl we provide the math ematical details of the transition layer (TL) 



model (ITitarchuk. Lapidus fc Muslimov Ill998l . hereafter TLM98) which is the basis for the 



proposed scenario. We present theoretical arguments for the inverse proportionality of QPO 
frequency on BH mass in §5.1[ We discuss the signature of bulk motion Comptonization 
(BMC) and its relation with the index evolution during state transition in §5.21 Specifically, 
we show that the index saturation effect is a direct consequence of the existence of the in- 
nermost bulk motion region and, therefore, can be considered as an observational signature 
of a BH. In ^we summarize our results and compare BH masses and distances with those 
obtained by other methods. We also discuss briefly the alternative methods for BH mass 
determination from X-ray data. Conclusions follow in [|7l 



2. Observations and data reduction 

For the study presented in this Paper we have analyzed data for seventeen spectral 
transitions in seven BH transient sources observed with RXTE. We also analyzed a subset 
of Cygnus X-1 data taken during 2001-2002 years when the source showed a transition from 
LHS to HSS. In Tabled] we summarize data for each transition and also provide references 
to previous analyses. 

We use the archival RXTE data from the HEASARC0. We present RXTE/PCA light- 
curves and hardness ratios for all analyzed outbursts in the top panels of Figures [2ll6l These 
diagrams are created using RXTE mission long data products provided by RXTE team. 
Each data point is based on a separate pointed PGA observation. Each transition data set 
is indicated by a separate color. Black color indicates observations during HSS or extreme 
LHS which are not used in our analysis. 



^http:/ /heasarc. gsfc.nasa.gov/ 
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RXTE/PCA spectra have been extracted and analyzed in the 3.0-50 keV energy range 
using the most recent release of PCA response calibration (ftool pcarmf vll.l). We cal- 
culated the root-mean-square normalized Fourier Power Density Spectra (PDS) from high 
time-resolution PCA mode data in the frequency range 0.01-1024 Hz. The relevant dead- 
time corrections to energy spec tra and timing spectra have been applie d following "The 
RXTE Cook Book" recipes and iRevnivtsev. Gilfanov fc Churazov I (120001 ) correspondingly. 
The PCA energy spectra were modeled using XSPEC astrophysical fitting software. Spec- 
tral analysis was done using approach similar to that adopted in ST07 for Cygnus X-1 and 
GRO J1655-40 data. To fit PCA spectra we used the sum of the bm c component (Generic 
Comptonization model, see iTitarchuk. Mastichiadis fc Kylafis 1 119971 ) and a Gaussian with 
the energy ~ 6.5 keV, which is presumably due to iron emission line. This model was also 
modified by the interstellar absorption, using the wabs model in XSP EC, with a hydrogen 



column value fixed at the Galactic value given by nH HEASARC tool (jPickey fc Lockman 



1990l ) for a particular source and a high energy cut-off (highecut). The width of the Gaussian 
was constrained in the range 0.8-1.2 keV. The high energy cutoff component accounts for 
exponential overturn of the spectrum due to electron recoil effect. Systematic error of 1.0% 
have been applied to the analyzed spectra. 

The bmc model describes the outgoing spectrum as a convolution of input "seed" black 
body spectrum having normalization Nbmc and color temperature kT with a Green function 
for Comptonization process. Similarly to the ordinary bbody XSPEC model, the normaliza- 
tion Nhmc is a ratio of source luminosity to square of the distance 



N, 



bmc 



1039erg/s/V d 



lOkpc 



:i^ 



The resulting model spectrum is also characterized by a parameter log{A) related to a 
Comptonized fraction / where / = A/{1 + A) and the Green's function spectral index 
a = r — 1, where P is photon index. 

There are two reasons for using the bmc model. First, bmc by the nature of the model is 



tion but also via dynamic or 


Dulk motion Comptonization (see 


Titarchuk. Mastichiadis & Kvlafis 


1997; 


Laurent & Titarchuk 


1999; 


Shaposhnikov & Titarchuk 1 


20061. for details). The sec- 



ond reason is that bmc calculates consistently the normalization of the "seed" photons Nbmc, 
presumably originated in the disk. The relation of N^mc and its proportionahty to the 
mass accretion rate in the disk is a direct consequence of the accretion disk theory (see e.g. 



Shakura fc Sunyaev Ill973l ). The adopted spectral model decribes well the most data sets 
used in our analysis. The value of reduced x^-statistic Xred — X^/^dof, where Ndof is the 
number of degrees of freedom for a fit, is less or around 1.0 for most of the observations. 
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For a small fraction (less than 3%) of spectra with high counting statistic the value of xled 
reaches 1.5. However, it never exceeds a rejection limit of 2.0. 

The PDS continuum shape in LHS and IS has usually a band-limited noise shape, which 
is well presented by an empirical model Px ~ (1.0 + (x/x*)^)"*"". The parameter is related 
to the break frequency and 2 x in is the slope of the PDS continuum after the break. We 
use Lorentzian shape to fit QPO features. We quote the Lorentzian centroid as a QPO 
frequency. 

Our approach to the analysis of a BH state transition is following. We identify a 
complete spectral transition as a period when the source state changes from LHS to HSS (or 
vice versa) showing gradual transition through IS. In this Paper we analyze the transitions 
for which the frequency of RXTE/PCA pointing observations allows to infer the index-OPQ 
correlation patterns and to perform a comparative study. In Table [T] we provide the list of 
the transitions (outbursts) which we utilized in our study. In this Table we also introduce 
an ID name of each transition which consists of abbreviated source name, dash, a transition 
direction identifier ("R" for the rise and "D" for the decay), and the two last digits of the 
year during which the event occured. In Table [2] we present the details of the model fits to 
the spectral and timing data for each RXTE observation in the spectral transition. 

3. Correlations between spectral and timing properties during state 

transitions 

3.1. Phenomenology of the correlation behavior 

In Figures [2ll6] (top panels) we show the spectral evolution of a given source throughout 
each analyzed outburst by presenting the PCA lightcurve and PGA hardness ratio. In the 
bottom panels from left to right we plot the correlations of photon index F versus QPO 
frequency, bmc model normalization Nbmc and Comptonized fraction / = A/{A + 1). Data 
for a particular transition are distinguished using the same color legend for the correlation 
plots and the lightcurve diagrams. 

The general behavior of the correlations is following. The QPO frequency and normal- 
ization are positively correlated with index. Comptonization fraction and index are usually 
anti-correlated (see right bottom panels in Figs. [2ll6]). In the upper part of the correlation 
patterns we observe a saturation of index which is more apparent in the index-normalization 
domain than in the index-QPO pattern. The reason for this is that the index and nor- 
malization values, as the best-fit parameters of the energy spectrum, are always available 
independently of the source state, while QPOs tend to faint as a source enters HSS. 
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In ST07 we claimed the presence of index saturation in Cygnus X-1, GRO J1655-40 and 
GRS 1915+105 using correlations with QPO frequency. Here we confirm the index saturation 
effect for GRO J1655-40 and Cygnus X-1 (see Fig. [2]andE]). Moreover, we establish the 
index saturation effect in GX 339-4, H 1743-322, 4U 1543-47 and XTE J1650-500 (see Figs. 
|3]l5] and [8] respectively). It is important to mention that the saturation of index with respect 
to the photon spectrur n normalization (mass a ccretion rate) in XTE J1650-500 was also 



recently established by iMontanari et al. I (120081 ) in their analysis of BeppoSAX data from 
2001 outburst. 

Three rise transitions, J1550-R98, H1743-R03 and J1859-R99, exhibit similar behavior. 
For all these transitions, when the index reaches its high values (>2.5), the Comptonization 
fraction shows non-monotonic behavior. Namely, for J1550-R98 and H1743-R03, as the 
source leaves LHS / starts to decrease as expected. However, when the source enters the IS 
the index-/ track reverses and / returns back to its maximum value of 1.0. In Figure H] we 
illustrate this evolution of the Comptonization fraction / as a function of photon index F 
for the rise transition H1743-R03 (see red data points on the right bottom panel). 

During the decay transitions the source luminosity is lower than that for the rise episodes 
and the corresponding variabilit y is suppressed, which is probably due to lower excitation 



power in the accretion flow (see iTitarchuk fc Shaposhnikov I l2008l . for details). As a result 



the spectral and timing data related to the decay transitions have lower signal-to-noise ratio 
than those for the rise transitions. The index saturation levels for the decay stages have 
values close to 2.1. Most probably, this is due to the fact that at low accretion rates during 
an outburst decay the cooling of X-ray emission area [Comptonization cloud (CC)] by the 
disk photons is less efficient than that at high M. As a result, the CC plasma temperature 
Tg is higher, leading to more efficient Comptonization and harder spectrum. This direct 
connection between the temperature of t he CC and the hardnes s of th e emergent spectrum 



was shown by Monte-Carlo simulations in iLaurent fc Titarchuk I (119991 ). In other words, the 



decay transitions appear to be different from the rise transitions: correlation patterns are 
not scalable in the normalization domain and they are usually not self-similar given that the 
index saturations take place at different levels. Our experience shows that while decay-to- 
decay and rise-to-rise scaling produces consistent results, decay-to-rise pattern scaling is not 
a reliable mass-distance indicator. 

To illustrate the above points, let us consider properties of state transitions in GX 339-4. 
In this source five different (rise and decay) transitions exhibited five different tracks (see 
Figure [3], third and fourth row panels from the top). Two rise transitions, GX339-R02 and 
GX339-R07, occurred at higher fluxes than that for two decay transitions, GX339-D03 and 
GX339-D07. These rise transitions also showed steeper index-QPO (index-normalization) 
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correlations with higher saturation levels. The difference between these transitions is clearly 
seen. On the other hand, the GX339-R04 hard-to-soft transition resembles a decay track 
with the index saturated at the value of 2.15 and a lower luminosity. But even for this rise 
episode the index-normalization track is much steeper than those for the decay episodes. The 
decay transition during 2004 outburst was quick and did not show sufficient QPO activity. 
We, therefore, discarded this transition episode from our analysis. 

It is also worthwhile to discuss the behavior of XTE J1550-564 during the 1998 rise 
outburst in detail. The outburst started on MJD 51063 and went through the initial LHS 
phase, which ended on MJD 51069. During this period spectral index increased from 1.35 to 
2.2 and QPO frequency changed from 0.3 to 3.3 Hz. In Fig. [7] this stage is marked by filled 
black circles. For the following 10 days the source exhibited IS with a strong fiare on MJD 
51076, when the index and QPO frequency peaked at 2.8 and 14 Hz correspondingly (IS-VHS 
track). We mark these data points by filled red circles. After the fiare the source returned 
to the IS with index ~ 2.3 and QPO at ~ 3 Hz. For the next 20 days we observed smooth 
evolution toward the HSS with index slightly increasing to 2.5 and QPO frequency rising 
to ~ 5 Hz. On about MJD 51105 the source entered the HSS (IS-HSS track). The index- 
QPO correlation consists of two different tracks which are shown on Fig. [3 The correlation 
curve related to the IS-HSS stage is lays lower than that related to the IS-VHS transition. 
Clearly, the conditions in the accretion flow changed between these episodes. During the 
IS-VHS stage the source presumably went through the strongest surge of accretion. This 
cold accretion flow provided strong cooling, which is manifested by higher spectral index. 
These changes in the physical conditions in the accretion flow led to a double-track pattern 
of the correlation. Therefore, for the scaling purposes we treat a combination of the initial 
LHS with IS-VHS and IS-HSS data as two separate correlation patterns. We further refer 
to them as J1550-98Ra and J1550-98b respectively. 

We revisited the data for GRO J1655-40 and Cygnus X-1, which were used for BH mass 
measurement in Cygnus X-1 (ST07). We present the data for these sources in Figs. [2] and 
ini In fact, Cygnus X-1 is a persistent source, with somewhat different phenomenology from 
the other sources analyzed in this study. The Cyg X-1 data obtained during 2001 and 2002 
were chosen for the analysis (see Fig. [H]). In the 2001-2002 period the source exhibited two 
transition, from LHS to HSS (red data points) and a reverse transition marked by blue color. 
Black points indicate the observation when the source was in the HSS. As seen from the lower 
panels, the rise and the decay data from Cygnus X-1 are very similar in terms of the index- 
QPO and index-normalization correlations. In particular, the same index-normalization 
tracks indicate that there is no hysteresis effect in this source. This may be related to the 
persistent nature of the accretion or the 0-star wind. We, thus, use both the rise and the 
decay data for Cyg X-1 data as one data set to construct the correlation pattern. This 
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approach is consistent with that adopted in the ST07 analysis. 



3.2. The scaling technique and parametrization of the correlation patterns 

One of the goals of the presented study is to apply X-ray observational data from BH 
sources to infer their fundamental characteristics. To do this we use the BH state paradigm 
combined with the scaling laws expected to be observed in spectral and timing data during 
state transitions. 

First, using available X-ray observations we verify a theoretical conjecture that the 
QPO low frequency is inversely proportional to a BH mass m = Mbh/MQ (see TLM98 and 
TF04). This QPO frequency dependence on m is a simple consequence of the fact that all 
characteristic dynamical tim escales in an accreting flow onto a compact object are set by the 



mass of the object (see e.g. IShakura fc Sunyaev Ill973l . and §5.11) . Thus, QPO frequencies 



being in the inverse proportion to these timescales, have to be inversely proportional to m. 
This intuitive law of nature now finds its observational confirmation. In §5.11 we present 
theoretical arguments for the inverse dependence of QPO frequencies on m in terms of 
transition layer model (see details in TLM98, TF04 and Appendix A). 

Thus, we write the first scaling law (that we implement in our scaling technique for the 
BH mass determination) in terms of a scaling factor 

s ='± = ^ (2) 

where subscripts r and t of frequency u denote reference and target sources respectively. 

The second scaling law, which we use as a basis for our mass and distance determination 
analysis, relates a source fiux F detected by an observer on the Earth, a disk fiux L and 
distance d, namely: 



The disk fiux L can be represented as 



(3) 



r]{r^)mLEd- (4) 



Here i?* = r^,i?g is an effective radius where the main energy release takes place in the disk, 
Rs = 2GM/c^ is the Schwarzschild radius, t] = l/(2r*), rh = M/ Merit is dimensionless mass 
accretion rate in units of the critical mass accretion rate Merit = -Z^Ed/c^ and Leci is the 
Eddington luminosity. 
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(5) 



i.e. M and thus using Eqs. |31l5]we have that 



(6) 



Note that both m and t] are assumed to be the same for two different sources in the same 
spectral state, which leads to Lr/Lt = Mr/Mt = l/s^. In our spectral analysis we determine 
the normalization of "seed" radiation, which is supplied by an accretion flow (disk) prior to 
Comptonization. The ratio of these normalizations for target and reference sources in the 
same spectral state can be written as 



Here fa is a geometrical factor which takes into account the difference in viewing angles of 
the seed photon emission area (the inner disk) by the Earth observer for two given analyzed 
sources. Therefore, in the case of radiation coming directly from the disk it would have the 
value /g = (cos^)r./(cos^)t, where 6 is the inclination angle of the inner disk. Note that 
this form of fc, however, may not be very accurate because the geometry of both the inner 
disk and the corona can be different from the plane (disk) geometry. In fact, in BH states of 
interest (LHS ans IS) they can be described, for instance, by quasi-spherical configuration. 
Despite of this uncertainty in the determination of fc, we adopt the above form fc in which 
^ ~ z if the information on the system inclination angle i is available. In doing this we 
introduce an additional assumption that the inner accretion disk and the system inclination 
angles are the same. This assumption, however, proved to be in good agreement with our 
analysis. 

Now we write the final equations of our scaling analysis. Namely, when Si, and sn are 
measured, the mass and the distance of the target source can be calculated as 



The part of the scaling method presented by equation is new. In the previous analysis 
(see ST07) we only use the scaling factor of the index-QPO correlation s^, and apply it 
for determination of the BH mass distribution. The new modified technique uses timing 




(7) 



Mt = SyMr 




and 




(9) 
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information and the normalization of photon spectrum to constrain both BH mass and the 
distance for a target source, provided that these quantities are known for a reference source. 

With equations and in hand, the task of BH mass and distance measurements for 
a target source is reduced to the determination of scahng coefficients s^, and s at with respect 
to the data for a reference source. The appropriate technique has aheady been implementaed 
in ST07. Specifically, after scalable state transition episodes are identified for two sources, 
the correlation pattern for a reference transition is parameterized in terms of the analytical 
function 



where argument x is either QPO frequency u or normalization Nhmc- By fitting this functional 
form to the correlation pattern, we find a set of parameters A, B, D, xtr and (3 which represent 
a best fit form of the function / for a particular correlation curve. 

For X ^ Xtr the correlation function f{x) converges to constant value A. Thus, A is a 
value of the index saturation level, /3 is the power-law index of the part of the curve for lower 
argument values and Xtr is the frequency /normalization at which index levels off. Parameter 
D controls how fast the transition occurs. We scale it to a target pattern by a transform 
X ^ sx when the best fit is found for a reference pattern. 

The use of analytical function to parametrize the correlations allows us to avoid theo- 
retical modeling of the index-QPO relation based on any particular physical model, which 
would require some internal parametrization in any case (see discussion in ST07). It makes 
the scaling procedure model independent, and it is justified because it satisfactory represents 
both the index-QPO and the index-normalization patterns. 



The Galactic microquasar GRO J1655-40 is well studied BH X-ray binary with the best 
measured mass and distance among known stellar mass BH sources. We therefore use this 
system as a primary reference source for our scaling. In addition, the ancillary information 
on the system inclination is available which makes distance determination more reliable. 
However, not all sources we analyse here have transitions directly scalable to GRO J1655-40 
data. Therefore, we have to utilize intermediate reference source to provide measurements for 
all sources. Specifically, only GX 339-4 and 4U 1543-47 have transitions which are scalable 
to the our primary reference source. After measuring BH mass and distance for GX 339-4 
we use this source as a reference for XTE J1550-564 and XTE J1650-500. Note, this scaling 
between GX 339-4 and XTE J1650-500 was allowed both with rise and decay data. We 




(10) 



4. 



BH mass and distance measurement in XRBs 
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then use XTE J1550-564 as a reference source for scaling to H 1743-322 and XTE J1859-226 
data. Finally, we measure BH mass and distance for Cygnus X-1 by scaling to GRO J1655-40 
and 4U 1543-47. All measurements are summarized in Table [3] and illustrated in Figures [H] 
and M Below we provide some details on how we combine our scalings in chains to provide 
measurements for all sources of interest. 

We first determine parameters for sources which have correlation patterns self-similar 
to GRO J1655 data. Namely, data sets GX399-D07 and U1543-D02 are scalable to that of 
J1655-D05. These scalings provide BH mass measurements for GX 339-4 and 4U 1543-47 
(Fig. E] panel A and B). Other decay transitions (i.e. GX339-D03, J1650-D01 and J1743- 
D03) have slightly higher saturation levels of ~ 2.1 and present a different cluster of scalable 
patterns (datasets). Thus we use GX339-4 as a reference source to provide the mass-distance 
determinations for XTE J1650-500 (Fig. [S] panel D) and H 1743-322 (Fig. [g] panel F). We 
use GX339-R04 data to scale to XTE J1550-564 (J1550-R00) data and to measure BH mass 
and distance in XTE J1550-564 (Fig. [8] panel C). The resulting value of BH mass is in an 
excellent agreement with the dynamical determination of BH mass in XTE J1550-564. We 
also scale the rise transition datasets GX339-R04 and J1650-R01 (Fig. [8] panel E) to provide 
a double check with the decay data scaling (see Fig. [S] panel D) for the same sources. The 
results obtained with these parallel scalings are in a striking agreement. 

Finally, we use XTE J1550-564 (J1550-R98) data for scaling with H 1743-322 (1743- 
R03) and XTE J1859-226 (J1859-R99) (Fig. panels G and H correspondingly). As we 
mentioned in the previous section, during the 1998 outburst rise from XTE J1550-564 the 
index exhibited two correlation tracks with different saturation levels. The initial LHS 
combined with subsequent IS-VHS phase is similar to H 1743-322 transition observed during 
the 2003 rise episode (1743-R03). We therefore use this subset of XTE J1550 (J1550-R98a) 
to scale to the 1743-R03 data. The IS-HSS phase of J1550-R98 data is scalable with J1859- 
R99 transition (Fig. [H panel H). We also note that these outbursts (J1550-R98, J1859-R99 
and H1743-R03) are similar not only in terms of the index-QPO and index-normalization 
behavior but also in terms of the index-/ patterns and the overall appearance of outburst 
lightcurves (we illustrate these observational appearances in Fig. H] which presents data for 
H 1743-322). 

Finally, in Fig. [9]panel K, we show the scaling between J1655-D05 and CYGXl datasets 
(see Tabled] for these dataset descriptions). In Table [3] we present the scaling coefficients Si, 
and Sn obtained using our scalings procedure. We also provide BH masses and distances 
obtained using the calculated scaling coefficients and Eqs. dHJED- The parameterizations of 
reference patterns in terms of the empirical functional form f{x) (Eq. [TOj) are presented in 
Table H 
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Thus, we construct the sequence of scahngs which allows us to obtain the BH mass and 
the source distance for a number of BH XRBs (and in some cases to provide a double check 
for these measurements). Table presents the final results of BH mass and distance de- 
terminations using our scaling method along with independently reported masses, distances 
and inclination angle values found in literature. For the sources for which we were able to 
make parallel scaling measurements with two different dataset pairs, the final best-fit values 
are the averages between individual measurements and the error ranges are calculated using 
the sum of error intervals given by individual scalings. 



5. Physical scenario for the evolution of spectral and timing properties during 

BH state transitions 

The accretion process onto a BH is far from being completely understood. Theoreti- 
cal models which are available for astronomers to explain observational phenomena usually 
deal with only narrow aspects of a broad observational picture. For example, while t here 



are several proposed QPO models in the literature (see iRemillard &: McClintockl |2006| . for 
references), most of them concern only QPO phenomenology and lack the connection with 
spectral properties and state transition aspects. However, in the recent years a concept of 
"transitional layer" (TL) in the accretion flow has emerged. In the framework of this model 
several major observational aspects of accreting BHs find natural explanation. We argue that 
this concept is a strong candidate to provide a basis for a theory of accretion process onto 
compact objects which would accont for the current observational picture of BH sources. 
Moreover, this paradigm may provide us with a first direct BH signature observationally 
confirmed from multiple sources. 

The starting point in the development of the TL paradigm is the notion of the necessary 
deviation of the rotational profile in the innermost part of the disk set by the conditions at (or 
near) a compact object for which rotation is presumed to be much slower than the Keplerian 
frequency near the object. This problem was, the first time, considered by TLM98 where the 
authors showed that the adjustment of the Keplerian rotation in the accretion disk to the 
sub-Keplerian rotation of the central object leads to the formation of the inner hot Compton 
corona (CC). 
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5.1. The inverse dependence of the QPO frequency vl on the BH mass 

TLM98 define the transition layer as a region confined between the the inner sub- 
Keplerian disk boundary (where the inner TL radius Rin = bRs, Rs is the Schwarzchild 
radius and is a numerical factor ^3) and the first Keplerian orbit (see Figs. 1 in TLM98 
and TF04). Numerical calculations (see also Eq. IA4I and Fig. 2 in TLM98) showed that 
the value of dimensionless outer radius Vout = Rout/ Rin strongly depends on 7— parameter 
(Reynolds-number) when the rotational velocity of central of object uo is less then the 
Keplerian frequency cjk at the TL outer boundary. 

The fact that Tout is a function of 7— parameter implies that the CC dimensionless size 
Icc = {Rout — Rin) /Rin = rout — 1 is a fuuctiou of 7 too. Specifically, given that 7— parameter 
defines the spectral state (see TLM98 and TL04) the CC dimensionless size Icc is the same 
for any given spectral state of BH even if BH masses differ by orders of magnitude. Thus the 
CC dimensional size Lcc = bRslcdl) ~ 9/cc(7)^ km is proportional to BH mass m = M/Mq 
for a given spectral state (or 7) . Moreover TLM98 demonstrated that the dimensionless CC 
size lcc = Lcc/ (bRs) anti-correlates with the X-ray spectral hardness (spectral index). TF04 
identified the QPO frequency ul as a normal mode oscillation frequency of the CC near the 
BH. Thus one should expect that the QPO frequency is inversely proportional to the BH 
mass given that the QPO frequency is proportional to a ratio of the characteristic velocity 
with which the perturbation propagates in the CC, note it is a function of spectral index 
only, and the CC size Lcc oc hcRs- The CC plasma (sonic or magneto-acoustic) velocity Vma 
decreases with M, related to an increase of 7, because in the case of the high mass accretion 
rate the disk provides strong soft photon illumination of the CC area. In other words, Vma 
is a strong function of the spectral state (or 7). Given that the QPO frequency should 
be expressed as a ratio of Vma and the CC size Lcc, namely 

oc - (11) 

Lcc{l,rn) dmlcciv km 

one can conclude that (i) there should be a correlation of with spectral state (or 7) and 
(ii) for a given spectral state z/^ should be inversely proportional to BH mass m. Thus, if 
this index-QPO correlation occurring during a spectral transition has a form which is similar 
to that for another BH source then there is a possibility to determine a ratio of BH masses 
for these two BH sources by sliding their correlations along log axis with respect to each 
other (see ^for details of an implementation of this scaling technique). 

R ecent detailed analysis of the spect ral and timing properties of X-ray emission in Cyg 



X-1 by lTitarchuk fc Shaposhnikov I (120081 ) confirms that QPO frequency indeed correlates 
with the Reynolds number 7 and photon index F. The Comptonization Y parameter (or the 
photon index F) which characterizes the efficiency of the Comptonization (see more details 
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in Appendix B) is determined by the product of the average fractional energy change per 
scattering rj and mean number of scattering N^c in the CC. TF04 used the results of the 
Monte Carlo simulations of Comptonization processes in the CC and the bulk motion flow 
to demonstrate that F increases and then saturates when M increases. Thus given that 
and r correlate with M the QPO frequency vl should correlate with the photon index T as 
well. 



5.2. The index saturation vs disk flux as an observational evidence of the 

converging flow (bulk motion onto BH) 

As we have already pointed out, positive correlations between photon index V and QPO 
frequency should be observed during the state transition, when the corona is cooled down by 
the disk photons (TLM98, TF04). As the temperature of the CC drops during the transition 
towards softer states, the bulk motion Comptonization (BMC) effect becomes dominant in 
the formation of the hard tail of X-ray spectra in BHs in the final spectral stage, i.e in HSS 
(see 



Titarchuk fc Zannias Ill998l . hereafter TZ98, and LT99). 



As seen from Fig. [2l-l6] the index saturation is detected in index-6mc normalization 
correlations for most analyzed sources (sometimes the index saturation is seen in the index- 
QPO correlations too). Namely, the index saturation effect is found in GRO J1655, GX 
339, XTE J1550, XTE J1650, H 1743, 4U 1543, Cyg X-1 (see Figs [Ml and [8] respectively). 
The saturation level of index can vary from source to source. Even for the same source the 
index saturation can be variable for different transition episodes. For example GRO J1655- 
40 exhibits the different index saturation levels Vgat ~ 2.3, and 2.05 during 2005 outburst 
rise and decay stages respectively (ST07). 

This effect of index saturation vs optical depth of the bulk flow (BM) r was first predicted 
by TZ98 and subsequently reproduced in Monte-Carlo simulations by LT99. It is worth 
noting that the index saturation effect is an intrinsic property of the bulk motion onto BH 
given that spectral index a = F — 1 is a reciprocal of the Comptonization parameter Y (see 
Eq. IB7p which saturates when the BM optical depth r, or dimensionless mass accretion 
rate m, increases. In fact, Y-parameter is a product of the average fractional photon energy 
change per scattering rj and the mean number of photon scatterings Ngc-, i-e. Y = rjNsc- 

The thermal Comptonization parameter Y ^ (4fcT/mpC^)r^ given th at in this case 



rj = ikT/mec"^ and Ngc ~ for r 1 (see e.g. iRybicki &: LightmanI 119791 ) and, thus, the 
thermal Comptonization spectral index 

a ~ [{AkT/meC^)T'^]-\ (12) 
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In the case of converging flow the preferable direction for upscatterred photons is the direction 
of bulk motion onto the BH, i.e along its radius. Note the fractional photon energy change 
is 

^E/E = (1 - ^llVR/c)/{l - /is^c). 

where /xi and /i2 are cosines of angles between the direction of electron velocity n = ^ r/Vr 
and direction of incoming and outcoming (scattered) photons respectively. 

The number of scattering of up-Comptonized photons A^^c can be estimated as a ratio 
of the radial characteristic size of the converging flow L and the free path / in the direction 
of motion, namely Ngc oc L/l = t given that AE/E has a maximum at /i2 = 1 for given /ii 
and Vr . On the other hand the effi ciency per scattering for bulk motion flow t] oc l/r when 



r > 1 (ILaurent fc Titarchuk 1120071 ). 



Thus one can conclude that the Comptonization parameter Y = rjNgc and energy index 
a = (see Eq. IB7|1 saturate to a constant values when optical depth (or mass accretion 
rate) of the BM flow increases. 

The index saturation value is determined by the plasma temperature during a transition 
[see LT99]. On the other hand, (see e.g. TLM98) the plasma temperature strongly depends 
on the mass accretion rate in the bulk motion region Mhm and its illumination by the disk 
photons Fdisk- For higher M^m and Fdisk the plasma temperature is lower. The level of 
the index saturation decreases when the plasma temperature in the bulk motion increases 
(TF04). Thus the index saturation levels can be different from source to source depending 
on the strength of the disk in a particular case. 



6. Discussion 

Currently a BH identiflcation in X-ray observational astronomy is based solely on a 
mass of compact object. Namely, a compact X-ray source is classifled as a BH if it is 
well established that its ma ss exceeds the upper ma ss limit of a stable rotating neutron 



star, namely 3.2 (see e.g. iRhoades &: Ruffini II 19741 ). To the date there is only one widely 



accepted method for the mass determination related to the measurement of the mass function 
/(M) based on optical spectroscopy. While the theoretical mass function is a function of 
the two masses and the inclination angle, the value of the observationally inferred mass 
function (which is the minimum possible mass of the compact object) depends only on the 
velocity half-amplitude K and the orbital period P and independent of the inclination angle 
i : /(M) oc PK^. 
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Qrosa (120031 ) . has summarized the measurements of the rotational velocities and incli- 
nations for 17 black hole binaries. He comments that the difficult part of this procedure has 
been related to the inverse problem of parameter optimization since the parameter space to 
search is usually vast and often parameters can be tightly correlated. In many cases, the 
observed light curve is not entirely due to the companion star owing to the presence of extra 
sources of light (usually from an accretion disk). Another difficult aspect is the uncertainty 
in radial velocity amplitude, which appears as third power in solution. Therefore, while the 
dynamical mass measurements in a number of BH sources were successful, one should always 
be cautious about potential systematic errors. 

The dynamical measurement of the mass of a compact object of more than the theoret- 
ical limit for a stable NS configuration is currently considered a sufficient evidence that the 
object is a BH. We note, however, that by itself such a mass measurement is not a direct 
proof of the nature of an object as a BH. Such a proof, for example, would come from the 
confirmation of the presence of a BH horizon. Direct observations of a BH horizon, by theory 
are not possible. Therefore, one has to look for signatures which would require the BH hori- 
zon and therefore would serve as a proof of its presence. In this work we present theoretical 
arguments and observational evidence in support that the index saturation effect observed 
during BH spectral transitions is due to the dynamical Comptonization in the converging 
bulk motion in the innermost part of a compact object. Formation of this bulk motion region 
is not possible in the presence of any solid surface, which would otherwise be manifested by 
observations of a "feed-back" in the form of a strong spectral component due to the energy 
release on the surface and/or coherent pulsations as in case of neutron stars. In fact, none of 
this "feed-back" features were ever observed from these sources. Therefore, we argue, that 
the photon index saturation with the mass accretion rate is the signature of the bulk motion 
(converging) flow onto BH. In this Paper we present the observational evidence of the index 
saturation detected in several X-ray outbursts from BH sources. 

We measure the BH mass and distance using scaling relationships in the correlations 
between spectral and timing properties observed during the state transitions in Galactic 
accreting BHs. We complement previously proposed technique for BH mass determination 
using index-QPO shift (ST07) with scaling the index— X-ray spectrum normalization which 
contains information on a source distance. Using our modifled scaling approach we calculated 
BH mass and distances for seven sources. Results of these measurements are given in Table 
m We also present constraints on BH masses and distances which are available using other 
methods. 



First, we compare the BH masses found using dynamical methods with masses derived 
by the scaling technique. BH masses provided by scaling are in good agreement with those 
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given by the dynamical method for 4U 1543-47 and XTE J1550-564. For XTE J1859-226 and 
Cygnus X-1, BH mass constraints based on optical measurements are consistent with BH 
mass values obtained using the scaling technique. The BH ma ss of 12.3 =fc 1 .471^1^ obtaine d 
for GX 339-agreement with the lower limit of QMq provided by iMufioz-Darias et al.l (120081 ) . 
For H 1743-322 mass measurer nent by the dynarn i cal rn ethod is not available. Based on the 
values of high frequency QPOs lMcClintock et al.l (120071 ) concluded that the H 1743-322 has 
a high inclination {i > 70°) and the BH mass about 11 solar masses. Our BH mass values 
of m = 13.3 ± 3.2 agrees with this estimate. 

For XTE J1650-500 we obtain a new value for BH mass of m = 9.7 ± 1. 6 , whi ch is 
considerably higher than the the upper limit of 7.3 Mq inferred by lOrosz et al.l (120041 ). In 
this paper we als o retract our previous clairn of th e measurement of BH mass of 3.8 ± 0.5 in 
XTE J1650-500 (IShaposhnikov fc TitarchukI 120081 ). This preliminary claim for XTE J1650- 
500 was based on self-similarity of XTE J1650-500 and XTE J1550-564 index-QPO data 
only. In our study of index correlations with spectrum normalization and other parameters 
it has become clear that the data sets used to obtain low mass value for XTE J1650-500 
(J1650-R01 and J1550-98a) are not scalable. The new value of BH mass for XTE J1650-500 
is obtained by scaling with GX339-4 us i ng mo re careful consideration of scalability criteria. 
It is also worth noting that lOrosz et al.l (120041 ) remarked that their analysis were affected by 
poor statistic and limited number of templates. Moreover, their conclusion on the BH mass 
upper limit was based on the number of assumptions about the spectral class of the optical 
companion, accretion disk contribution and inclination angle. 

We hope that future observations will clarify the discrepancy in XTE J 1650-500 mass 
determinations by different methods. We, however, note an interesting fact that most our 
inferred BH masses, as well as majority of BH masses measured by dynamical methods fall 
in a narrow range around 10 Mq. This may be due to some observational selection effect or 
some currently unknown factors in stellar evolution. In our opinion, this matter deserves a 
separate investigation. 

In most cases our distance estimates are also in agreement with distances given by 
independent methods. Especially encouraging is a good agreement for 4U 1543-47 and 
Cygnus X-1. For this sources the inclinations are much smaller than that for the reference 
source GRO J1655-40. This means that an approximation of geometrical factor fg as a ratio 
of cosines of inclination angles works well. For H 1743-47 iMcClintock et al.l (120071 ) inferred 
the inclination angle of i ~ 70° and the distance of 10 kpc. Using this inferred value of i we 
obtained that fc ~ 1, and the distance of d = 9.1 ± 1.5 kpc for H 1743-322. 

An X-ray spectro scopic method of the BH mass determination was im plemented by 
Borozdin et al. I (119991 ) and further developed in IShrader fc Titarchuk I (119991 ). The method 
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is applicable to the HSS data when the thermal blackbody-like (BB) component, which is 
presumably related to the accretion disk emission, dominates in the emergent X-ray spec- 
trum. The main idea of this X-ray spectroscopic method is to use the color temperature and 
normalization of the BB component to infer the BH mass. One needs to employ a correct 
color factor Th which is a ratio of the color and effective BB temperatures of the disk. For the 
color factor calibration iBorozdin et al. I (119991 ) chose GRO J1655-40 for which d, m and i are 
well determined. Analyzing the RXTE syect ia. of this source and using the Shakura-Sunyaev 
(SS) disk model (IShakura fc Sunyaev Ill973l ) they obtaine d T/, ~ 2.6. This is drasticall y dif- 
ferent from the commonly used value of 1.7, obtained by IShimura &: Takahara I (119951 ) also 
in the framewo rk SS model but with certain assu mptions regarding the disk viscosity pa- 
rameter. Thus iBorozdin et al. I (119991 ) and later IShrader fc Titarchuk I (119991 ) used GRO 
J1655-40 calibrated value of the color ratio of 2.6 for the BH mass determination. The latest 
results of the BH mass determination f or a number of the BHC sou rces using the X-ray 
spectroscopy method are summarized in IShrader fc Titarchuk I (120031 ) . hereafter ST03. We 
also present these values in Table [5] for comparison. From this comparison one may see that 
the main uncertainties of the ST03 BH mass estimates, which relied on the spectroscopic 
method, are mostly driven by uncertainties of the inclination i, the distance d and by the 
value of the color factor Th. 



7. Conclusions 

We perform a comprehensive study of a highly representative set of well observed spec- 
tral transitions in Galactic BH sources. The goal of the study is to further explore the 
possibility of measuring the BH fundamental properties from X-ray data and search for BH 
signatures. We use correlations between spectral and timing properties during the state 
transitions as a main tool to investigate the BH spectral transition phenomenology. 

We compiled the state transition data from eight galactic BH sources collected with 
the RXTE mission. We examine the correlation between photon index of the Comptonized 
spectral component, its normalization and the QPO frequency. Analyzing the behavior 
of the correlation patterns we utilize four basic scaling laws: the inverse proportionality of 
frequencies of oscillatory processes to the mass of central BH and the disk flux proportionality 
to the BH mass, mass accretion rate and its inverse proportionality to a square of distance. 
We establish that scalable correlation patterns indeed contain information on BH mass and 
system distance. 

We then combine the scaling of the correlation patterns in the frequency and normaliza- 
tion domains for a set of Galactic X-ray BH binaries from which we determine BH masses and 
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distances. Our results confirm that the correlation scaling method is a powerful technique 
for a BH mass determination. The application of the scaling technique for the high precision 
measurements of BH masses requires very well sampled observations of the source evolution 
through the outbursts and careful consideration of scalability of correlation patterns. 

We have tested the scaling method using the known (from optical, IR observations and 
X-ray spectroscopic measurement) BH masses in 4U 1543-47, XTE J1859-I-226, XTE J1550- 
564 and Cygnus X-1. Using the inverse proportionality of QPO frequency to a BH mass 
for a given spectral state we arrive to a set of BH mass determinations, which are in good 
agreement with available dynamical data for several systems. 

We also show observationally that the transition layer theory (TLM98, TF04) correctly 
predicts the QPO frequency dependence on the BH mass and overall shape of the correlations 
between spectral and timing properties observed during spectral transitions. The success of 
the scaling method for the BH mass determination strongly supports the Compton Cloud 
(CC) origin of the observable index-QPO correlation. 

We present the observational evidence supporting the theory of the bulk motion (con- 
verging) fiow, i.e. the index saturation with the mass accretion rate. Specifically, analyzing 
RXTE observations for a set of BH outbursts we find that the index saturation is seen in 
the index-QPO and index-normalization correlations. We argue that the index saturation 
with mass accretion rate as a signature of the bulk (converging) fiow should only exist in 
the BH sources. Only in these sources there is no radiation pressure feedback effect at high 
mass accretion rate, as that takes place in NS binaries. In other words, this index saturation 
effect provides a robust observational evidence for the presence of black holes in these BHCs. 

The authors acknowledge the referee for many constructive suggestions to improve the 
paper quality and its presentation. 



A. The transition layer paradigm 



The radial motion in the disk is controlled by friction and the angular momentum 
exchange between adjacent layers, resulting in the loss of the initial angular momentum by 
accreting matter. The corresponding radial t ransport of the angular mom entum in a disk is 
described by the following equation [see e.g. IShakura fc Sunyaev I (119731 )]: 



(Al) 
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where M is the accretion rate in the disk and Wr^p is the component of a viscous stress tensor 
that is related to the gradient of the rotational frequency uj = 2ttv, namely, 

Wr^ = -27]HR^, (A2) 

where if is a half-thickness of a disk and rj is turbulent viscosity. The non-dimensional 
parameter that is essential for equation (lAip is the Reynolds number for the accretion flow, 

7 = M/An7]H = RVr/D, (A3) 

where K- is a characteristic velocity, and D is the diffusion coefficient. D can be defined as 
D = Vtlt/3 using the turbulent velocity and the related turbulent scale, respectively or as 
D = Dm = (? I (y for t he m agnetic case where a is the conductivity [e.g. see details of the 



D-definition in iLang I (119981 )]. 



Equations uj = ujq sX R = Rm = bRs (at the inner TL radius), a; = at i? = Rout 
[the radius where the transition layer adjusts to the Keplerian motion for which uk = 
{GM/R^f% and dw/dR = duK/dR at i? = Rout were assumed by TLM98 to be the 
boundary conditions. Note that here we set the inner boundary at Rin = bRs with b of 3. 
This value of Rin is valid for a BH whose spin a is less than 0.8. 

Thus, the profile u!{R) and the outer radius of the transition layer are uniquely deter- 
mined by the boundary conditions and the angular momentum equation (lAmA2|) for a given 
value of the Reynolds number 7 (see Eq. IA3p . 

The solution of angular momentum equation (]AltlA2p satisfying the above boundary 
conditions is equation (10) in TLM98 while the following equation (see also Eq. 11 in 
TLM98) 

3eout/2 = D,-fr^t + 2(1 - Di)r;„2 (A4) 

determines rout = Rout/ Rin as a function of 7— parameter. Here 9out = ^k{rout)/^o and 
Di = [9 out - rout)/{roZ - rout)- 

The adjustment of the Keplerian disk to the sub-Keplerian inner boundary creates 
conditions favorable for the formation of a hot plasma outflow at the outer boundary of the 
transition layer (TLM98), because the Keplerian motion (if it is followed by sub-Keplerian 
motion) must pass through a super-Keplerian centrifugal barrier region. 

Note one can see that the right hand side of Eq. (lA4p is linearly proportional to Bout 
if 9 out ^ 1 because Di oc 9out and \Di/rl^^\ ^ 1 for 9out ^ 1 and 7 > 2. As a result 
9out can be canceled from left and right hand sides of Eq. flA4p when 9out ^ 1- Namely 
the value of dimensionless outer radius Tout strongly depends on 7— parameter (i?e— number) 
and independent of 9 out for 9 out >> 1. 
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It implies that the CC dimensionless size Icc = {Rout — Rin)/Rin = rout — 1 is a function 
7— parameter (i?e— number) only if 6out ^1- It is worth noting that in the general case the 
size Icc is a function of 7 and Uq (or BH spin a) too. The direct scaling of Lcc = bRslcdl, 
with BH mass m is not possible anymore. There is a systematic shift of the values of Icdl, o,) 
for a given 7 because of BH spin a. 



B. Spectral index a as a reciprocal of the Comptonization parameter Y 

The intensity of the injected soft photons undergoing k scatterings in the Compton 
cloud is 

h^p' (Bl) 

where p is a mean probability of photon scattering in the CC. The probability of photon 
scattering is directly related to the mean number of scatterings 

00 

iVsc = 5ZA;/g = p/(l-p) (B2) 
fc=i 

where q = 1 — p is a probability of the photon escape from the CC. Thus, using Eq. (\B2\i . 
we obtain 

p = l-l/(iV,, + l). (B3) 

Because the average photon energy change per scattering < AE >= rjE (where t] > for 
upscattering case), the injected photon energy after k scatterings E is 

E = {l + r])''Eo (B4) 

The combination of Eqs. (IBip . (lB4p yields that the emergent upscattering spectrum of the 
soft photon of energy Eq in the bounded Compton cloud is a power law 

/.oc(|)"° (B5) 

which energy index is 

\n{l + ri) ^ ^ 

Using Eq. (1B3P we can reduce Eq. (1B6P to 

a ^ ivNscy' = Y-\ (B7) 

for Nsc > 1 and 77 < 1. 
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Table 1. Spectral Transitions Data Used in Analysis 



Source 


Transition ID 


Start 


End 


Type 


Refs. 


GRO J1655-40 


J1655-R05 


22/02/05 


13/03/05 


RISE 


11 




J1655-D05 


11/09/05 


27/09/05 


DECAY 


11 


GX 339-4 


GX339-R02 


23/04/02 


18/05/02 


RISE 


4,5 




GX339-D03 


20/02/03 


05/06/03 


DECAY 






GX339-R04 


10/07/04 


20/08/04 


RISE 






GX339-R07 


22/01/07 


18/02/07 


RISE 






GX339-D07 


07/05/07 


01/06/07 


DECAY 




4U 1543-47 


U1543-D02 


14/07/02 


31/07/02 


DECAY 


9,10 


XTE J1550-564 


J1550-R98 


08/09/98 


16/10/98 


RISE 


1,2,3 




J1550-R00 


10/04/00 


05/06/00 


RISE 


1,2,3 




J1550-D00 


05/06/00 


16/07/00 


DECAY 


1,2,3 


XTE J1650-500 


J1650-R01 


06/09/01 


28/09/01 


RISE 


7 




J1650-D01 


19/11/01 


25/11/01 


DECAY 


7 


H 1743-322 


H1743-R03 


28/03/03 


27/04/03 


RISE 


3 




H1743-D03 


15/10/03 


03/11/03 


DECAY 




XTE J1859+226 


J1859-R99 


09/10/99 


27/10/99 


RISE 


6 


Cygnus X-1 


CYGXl 


15/06/01 


01/12/02 


MIXED 
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Table 2. Spectral and Timing Characteristics for RXTE data^ 



TVansition ID Observation ID 


mjd"^ 


BH State 


^QPO 


a 


kT 


log{A) 


Nbmc 


Ecut 


Efold 


Fx X 10* 


2 

^red 




day 




Hz 




keV 






kcV 


keV 


erg/s/cm 




riv^^Q Rfi? Q9zi9R ni 09 nn 


r\AA 99 97 

04: IZZ . Z I 


LHS 




n A^^{Q\ 
u.^ioo^^y J 


n fii i''^'! 

u . u O j 


u.oy ^^o J 


n n'^fifi/'s^ 

U.UOUU^^O J 


Z U O ) 




n ^^79^9"! 
U.O 1 z J 


R9 
U.oz 




54127 11 


LHS 






U.OOl^O ) 


U.OOI^U ) 


n n7n^^9^ 

U.U / Ul^Z J 


Z O \^ I) ) 




1 .uoy i^y ) 


u.yo 


Q9n'^f^ ni ni nn 


1 OQ AT 
u'iJ.zy .'i ( 


T Nc: 


\(\A(f\\ 


n ^A {A \ 


n fi'^f''^^ 

U.Uo^^O J 


AAfO^ 


U.U ( U\^-L J 


21(2) 




1 ifin/'9^ 


n 7"^ 

u. / o 


090*^^^ m m m 

yzuoo-ui-ui-ui 


o^-Lzo.y^ 


T TT9 


149*'^^ 

U. L^£i\0 J 


A9i,A(9i\ 




1 Q('\'\ 
u. J-yi^o ) 


ORon 

U.UoU(^l J 


9^^ 






S7 

U.O 1 


Q9n'^ ^-01 -01 -09 


^4mi 10 


LHS 


lfifl(A'\ 
U. XOUl ^ ( 


u. oz 1 1 y^ 


u. oy \^o j 


U . OO ^^O J 


n C\R'^d(Q'\ 

U .yJOO'-tyiJ J 


26(1) 




1.276(3) 


1.03 


090'^'=; oi oi o'^ 


Cyl 1 Qfl 1 




U. ID ( \^ J 




u.ooi^z ) 


U.ZUI^O ) 




99/'i ^ 




1 914^''^^ 


U.OO 


090*^^^ ni ni (\A 


04:10Z.UO 


LHS 


9nnf''^^ 

U.ZUU^^O J 


U.OOl^ 1 J 


u.oyi^o ) 


u.oy ^^o ) 




24(1) 


1 1(0) 


1 314(2) 


1 n*^ 
1 .uo 


Q9n'^^-ni n9 nn 


•^41 '^'^ on 


LHS 


u.zo 1 yo J 


0.55(1) 


0.58(3) 


0.39(3) 


U . UOO ) 


25(4) 
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1.354(3) 


0.90 
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Zl .U^^O J 
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U. 1 z 


090*^^^ 01 0"^ 01 


f^AIAI 0^ 

O^J-^J. .Uu 


TS 


1 RQ1 f'i^ 


Q9n ^ 
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1 

U. lO^^Z J 


n 1 "ii fo^ 

U. lOl^Z ^ 


zu.u^^y J 






74 

U. 1 


92035-01-03-02 


54142.04 


IS 


2.434(9) 


1.03(1) 


0.51(1) 


-0.14(2) 


0.182(5) 


20.6(8) 


72(5) 


1.343(2) 


0.85 


92035-01-03-03 


54143.02 


IS 


3.508(9) 


1.22(1) 


0.55(1) 


-0.17(1) 


0.208(4) 


24(1) 


67(10) 


1.338(2) 


0.91 


92428-01-04-00 


54143.87 


IS 


4.34(1) 


1.29(1) 


0.597(9) 


-0.25(1) 


0.226(3) 


22(1) 


80(11) 


1.337(3) 


1.00 


92428-01-04-01 


54143.95 


IS 


4.24(1) 


1.26(1) 


0.52(1) 


-0.35(1) 


0.265(7) 


22(1) 


87(14) 


1.336(3) 


1.08 


92428-01-04-02 


54144.09 


IS 


4.13(2) 


1.29(1) 


0.57(1) 


-0.23(1) 


0.230(5) 


22(3) 


125(46) 


1.361(3) 


0.97 


92428-01-04^03 


54144.87 


IS 


5.00(2) 


1.35(1) 


0.59(1) 


-0.33(1) 


0.253(5) 


25(3) 


81(33) 


1.346(4) 


1.31 


92428-01-04^04 


54145.96 


IS 


5.63(2) 


1.39(1) 


0.597(9) 


-0.38(1) 


0.270(5) 


22(2) 


89(33) 


1.346(3) 


1.17 


92035-01-03-05 


54145.11 


IS 


5.80(2) 


1.35(1) 


0.554(8) 


-0.50(1) 


0.307(6) 


25(2) 


81(17) 


1.303(3) 


1.37 


92035-01-03-06 


54146.03 


IS 


5.55(2) 


1.36(1) 


0.601(8) 


-0.40(1) 


0.269(4) 


21(2) 


111(24) 


1.343(2) 


1.29 


92035-01-04-00 


54147.01 


IS 


6.67(2) 


1.53(1) 


0.626(7) 


-0.608(9) 


0.338(5) 






1.371(3) 


1.59 


92035-01-04-01 


54148.14 


IS 




1.26(3) 


0.796(7) 


-0.96(1) 


0.194(5) 






1.024(6) 


1.07 


92035-01-04-02 


54149.69 


HSS 




1.19(2) 


0.782(4) 


-1.104(8) 


0.193(4) 






0.963(6) 


1.27 



''Table in the paper/PDF version includes only data for GX339-R07 transition. Entire table is available in the electronic edition of the Paper. 
''Date at the start of the RXTE observation 
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Table 3. Index-QPO and Index- Norm Scaling Coefficients 



Target Tr. ID 


Reference Tr. ID 


Fig.E] 

Panel 


5. 


Sn 


m 


d 


fa 


GX339-D07 


J1655-D05 


A 


1.96±0.12 


1.65±0.06 


12.3±1.4 


5.75±0.64 


1 


GX339-R04 


J1655-R05 




1.95±0.19 


1.9±0.1 


12.3±1.8 


6.15±0.83 


1 


U1543-D03 


J1655-D05 


B 


1.49±0.15 


2.23±0.22 


9.4±1.4 


3.5±0.5 


0.37 


J1550-R00 


GX339-R04 


C 


0.87±0.02 


0.37±0.02 


10.7±1.5 


3.3±0.5 


1 


J1650-D01 


GX339-D03 


D 


0.80±0.02 


0.47±0.01 


9.9±1.4 


3.5±0.5 


1 


J1650-R01 


GX339-R04 


E 


0.76±0.02 


0.358±0.004 


9.3±1.2 


2.95±0.4 


1 


H1743-D03 


GX339-D03 


F 


1.15±0.07 


2.16±0.13 


12.7±2.6 


9.1±2.0 


1 


H1743-R03 


J1550-R98a 


G 


1.19±0.08 


6.44±0.52 


14.2±2.4 


9.1±1.6 


1 


J1859-R99 


J1550-R98b 


H 


0.696±0.015 


0.429±0.005 


7.7±1.2 


4.2±0.5 


1 


CYGXl 


J1655-D05 


K 


1.17±0.03 


0.144±0.003 


7.4±0.6 


2.04±0.18 


0.42 


CYGXl 


U1543-D03 




0.85±0.02 


0.064±0.001 


8.0±0.9 


2.1±0.4 


0.42 



''The column contains estimates for Jq based on available information on inclination angles. 



Table 4. Par ametrizat ions for reference patterns 



Transition 




Index-QPO 






Index-Normalization 








B 




D 




A 


B 


Ntr 


D 


/? 


J1655-R05 


0.59±0.04 


3.0±0.1 


1.0 


1.6± 0.2 


2.2 


0.400±0.003 


0.061±0.003 


1.0 


2.3±0.6 


J1655-D05 


0.58±0.01 


10.9±0.4 


0.1 


1.0 


2.02± 0.02 


0.44±0.02 


0.026±0.002 


1.0 


1.88±0.25 


GX339-D03 


0.66±0.02 


5.4±0.4 


0.1 


1.0 


2.08±0.01 


0.45±0.02 


O.OlOiO.OOl 


1.0 


2.0±0.3 


GX339-R04 


0.58±0.02 


1.4±0.2 


1.0 


1.0 


2.14±0.01 


0.51±0.02 


0.039±0.002 


1.0 


3.5 


J1550-R98a 


1.27±0.02 


1.84±0.07 


1.0 


0.65±0.02 


2.94±0.08 


1.8±0.3 


0.162±0.05 


1.0 


0.6±0.1 


J1550-R98b 


1.32±0.03 


1.97±0.07 


1.0 


0.5 


2.55 


1.25±0.04 


0.164±0.004 


1.0 


1.0 
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Table 5. BH Masses and Distances 



Source 










Mscal 


^scal ) 


Refs. 




Me 


deg 


kpc 


Me 


Mq 


kpc 





GRO n655-m'^ 


6.3±0.3 


70±1 


3.2±0.2 


GX 339-4 


> 6 






7.5±1.6 


4U 1543-47 


9.4±1.0 


20.7± 1.5 


7.5±1.0 


XTE J1550-564 


9.5±1.1 


72±5 


-2.5,~e 


XTE J1650-500 


2.7- 


^ 7.3 


> 50 


2.6±0.7 


H 1743-322 


~11 




~ 70 


~10 


XTE J1859-226 


7.6 - 


12.0 




11 


Cygnus X-1 


6.8 - 


13.3 


35± 5 


2.5±0.3 



1,2 



9.0±3.8{4) 


12.3±1.4 


5.75±0.8 


5,6 


14.8±1.6 (6.4) 


9.4±1.4 


9.4±1.8 


7,8 


9.4±2.1 (5) 


10.7±1.5 


3.3±0.5 


9,10,11 


10.6±4.0 (5) 


9.7±1.6 


3.3±0.7 


12,13 




13.3±3.2 


9.1±1.5 


14 


12.3±1.7 (11) 


7.7±1.3 


4.2±0.5 


15,16 




7.9±1.0 


2.2±0.3 


3,4 



^Dynamically determined BH mass and system inclination 
''Source distances found in literature 

'^Mi = Mcosi. Source distances used by ST03 are given in paranthesis. 

"^GRO J1655-40 is a primary reference source. All masses and distances given in columns 9 and 10 are 
determined with respect to the best measured parameters for this source. 
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Fig. 1. — Three representative energy (top) and power {bottom) spectra during the rise part 
of the 1998 outburst of XTE J1550-564. Data are taken from RXTE observations 30188- 
06-01-00 [red), 30188-06-01-00 {blue) and 30191-01-05-00 {black). In the top panel the sohd 
curves correspond to the best-fit model spectra. 
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Fig. 2. — Top: PCA Standardl count rate (left) and hardness (right) for the 2005 outbursts 
from GRO J1655-40. Blue and red colors stand for rise and decay transitions correspondingly. 
Botom: Photon index is plotted versus QPO frequency (left), BMC normalization (middle) 
and Comptonized fraction (right) for the rise (red) and the decay (blue) transitions. 
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Fig. 3. — PCA Standardl count rate {top row) and hardness {second row) for three out- 
bursts from GX 339-4 on 2002 (left), 2004 (middle) and 2007 (right). Third and bottom 
rows: Photon index plotted versus QPO frequency (left), BMC normalization (middle) and 
Comptonized fraction (right) for transitions in GX 339-4 (each transition is indicated by 
different color). 



- 33 - 



< 

o 
a. 



10000 



1000 



100 



10 



1 

52700 



"1 1 1 1 r 



J I I I L 



< 
O 



52750 



52800 52850 
MJD, day 



3 
2.8 
2.6 
2.4 
2.2 

2 
1.8 
1.6 
1.4 



~i — I — I — I — a — I — r 



t 



RISE 2003 
DECAY 2003 



J I I I I I L 



52900 52950 53000 



3.2 
3 
2.8 
2.6 
2.4 
2.2 
2 
1.8 
1.6 
1.4 
1.2 



0.5 
0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 



1 1 1 

_ + 


1 1 
















- + 


t\ 




++ +++ 








+ 

■ti- 






t 




1 1 "%*H-FtH^ 44-+-!^+^ 


+ 

1 1 





52700 52750 52800 



52900 52950 53000 



~i 1 1 1 1 1 r 



t i 



J I I I I I L 



4 6 8 10 12 
QPO Frequency, Hz 



14 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
BMC Normalization, Uq/d^m 




0.2 0.4 0.6 0.8 
Comptonized fraction, f 



Fig. 4. — Top: PCA Standardl count rate (left) and hardness (right) for the 2003 outbursts 
from H1743-322. Bottom Photon index plotted versus QPO frequency (left), BMC normal- 
ization (middle) and Comptonized fraction (right) for the rise and the decay transitions. 
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Fig. 5. — Top: PCA Standardl count rate (left) and hardness (right) for the 2002 outburst 
from 4U 1543-47. Bottom: Photon index versus QPO (left), BMC normalization (middle) 
and Comptonized fraction (right) for the decay transition. 
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Fig. 6. — Top: PCA Standardl lighcurve (left) and hardness (right) for the Cygnus X-1 
during 2001-2002. Red points indicate soft-to-hard transition while blue points show hard- 
to-soft transition. Observations during the HS are shown in black color. Bottom: Photon 
index versus QPO (left), BMC normalization (middle) and Comptonized fraction (right) 
during 2001-2002 state transitions in Cygnus X-1. 



- 36 - 




Frequency, Hz 




BMC Normalization, Lgg/d 

Fig. 7. — Photon Index - QPO frequency and Index-BMC normalization correlations (upper 
and bottom panels respectively) during the rise of XTE J1550-564 1998 outburst. Index- 
QPO correlation exhibits two different tracks, which correspond to different index saturation 
levels (see text for details). 
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Fig. 8. — Scaling photon index vs QPO frequency and normalization correlations for BH 
sources. Normalization plot is placed below the corresponding QPO frequency panel. In 
these diagrams the target transition for scaling is shown in red color while reference transition 
is plotted in black. Each pair (from left to right and from top to bottom) corresponds to the 
row in Table [3] beginning with the second row. 



- 38 - 




QPO Frgquency, Hz QPO Frequency, Hz QPO Frequency, Hz 




BMC Normalizaiion, Lgg/d^io BMC Normalization, Lgg/d^^Q BMC Normalization, Lgg/d^^Q 



g. 9. — Index-QPO frequency and index-normalization scaling diagrams continued. 



